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Abstract 

Constitutive equations are derived for the viscoelastic response of rubbery poly- 
mers at finite strains. A polymer is thought of as a network of long chains connected 
to temporary junctions. At a random time, a chain detaches from a junction, which 
is treated at transition from its active state to the dangling state. A dangling 
chain at a random time captures a new junction in the vicinity of its free end and 
returns to its active state. Breakage and reformation of long chains are modeled 
as thermo-mechanically activated processes. Stress-strain relations for a rubbery 
polymer are developed using the laws of thermodynamics. Adjustable parameters 
in the model are found by fitting observations in uniaxial tensile tests for a car- 
bon black filled rubber at various temperatures. Fair agreement is demonstrated 
between experimental data and results of numerical simulation. 

1 Introduction 

This paper is concerned with the nonlinear viscoelastic behavior of rubbery polymers at 
finite strains. The viscoelastic and viscoplastic response of rubbery polymers has been 
the focus of attention in the past decade p|-p9[|. This may be explained by the neces- 
sity to develop constitutive equations for rubber-like materials which may be employed 
for the numerical solution of applied engineering problems, on the one hand, and by a 
number of peculiarities in the behavior of rubber which cannot be adequately predicted 
by conventional stress-strain relations, on the other. The time-dependent response of 
industrial rubbery polymers in conventional mechanical tests is affected by a number of 
accompanying phenomena: 

1. Plastic deformation and micro-damage (initiation and development of shear bands 
and crazes) i, P, ||. 

2. Meso-scale damage (cavitation) and void formation |31| , [32]. 

3. Mechanically induced crystallization of semicrystalline polymers pl| , |2"7 , 3"3"| . 



1 



4. Slip |31] and fragmentation |35| of crystallites driven by network stretching. 

5. Mullins's effect g |, g. 

6. Payne's effect |T§. 

7. Strain induced evolution of the internal structure J/1] and orientational hardening 



21, 27, 36 



Reinforcement of rubbery polymers by carbon black and silica 0, H, [13], Dj[ [2^, [2S 
and fractal aggregation of particles JlO| . 



Although this list is far from being exhausted, even these items show that the existing 
micro-mechanical models are oversimplified to quantitatively predict the time-dependent 
response of rubbers. Even if we disregard these phenomena for an instant, a number of 
unresolved problems is found out which are associated with the description of rubbery 
polymers at the micro- level. 

This study focuses on the following questions: 



1. According to the entropic theory of rubber [38], a polymer may be treated as 



a network of long chains connected to junctions. A junction is thought of as a 
crosslink between two (or more) monomers belonging to different chains. A chain 
between two junctions consists of a large number of strands (statistically indepen- 
dent elements of a macromolecule) which are modeled as rigid rods. This implies 
that the mechanical energy of a chain vanishes and the main contribution to the free 
energy of a chain is provided by the configurational entropy (that characterizes the 
number of configurations available for a chain which is treated as a random walk on 
a three-dimensional lattice) ||38|| . This approach correctly predicts the response of 
polymeric gels near the sol-gel transition point. Its application to more dense media 
(e.g., to polymeric melts) results in large deviations between experimental data and 
results of numerical analysis. To avoid these discrepancies, a mechanical energy 
of chains is introduced "by hand" into the formula for the free energy (which 
is tantamount to treatment of strands as linear or nonlinear elastic springs). On 
the other hand, a conventional standpoint is that the concept of entropic elasticity 
is inapplicable to glassy polymers, whose response is associated with cooperative 
relaxation in micro-domains (regions of relatively high density comprising dozens of 
strands which belong to different chains and which rearrange simultaneously because 
of large-angle reorientation of strands) j|0[ [41], f|2|, [43| . An open question is about 



the correspondence between the entropic contribution into the free energy of a rub- 
bery polymer and the contribution of its mechanical energy at various temperatures 
T (in particular, near the glass transition temperature T g ). 

2. Changes in instantaneous elastic moduli of rubbery polymers with temperature are 
conventionally described by two mechanisms at the micro-level. According to the 
first J3£|, an increase in temperature results in thermal expansion of polymers, which 
is equivalent to an increase in the average distance between sequential monomers 
in a chain and between neighboring chains. The growth of this distance implies 
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a decrease in the intramolecular and intermolecular forces which is observed as a 
reduction of instantaneous elastic moduli at the macro-level. According to the other 
approach [fill , the rubbery state of a polymer is characterized by a large number 
of entanglements between long chains. At elevated temperatures, entanglements 
do not impose restrictions on micro- motion of individual strands (local constrains), 
but serve only as global topological constrains for the motion of macromolecules. 
With a decrease in temperature (or an increase in the strain rate), the number 
of entanglements grows that restrict the micro-motion of particular strands of long 
chains, which is equivalent to transformation of these entanglements into temporary 
crosslinks. An increase in the number of temporary crosslinks is tantamount to a 
decrease in the average number of strands in a chain. The question is which of these 
two mechanisms is more adequate for the description of observations in mechanical 
tests. 



3. According to conventional theories of rubber elasticity |37|, a polymer is treated as 



an ensemble of flexible chains (the energy associated with their bending is neglected). 
The account for mechanical energy of elongation of chains in the formula for the 
free energy implies the question whether a model of semiflexible chains (with finite 
bending rigidities) is more adequate for the description of rubbers, and in what way 
the bending rigidity should be introduced into the constitutive equations. 

This paper is an attempt to shed some light on these issues. To describe the time- 
dependent behavior of rubbery polymers, we employ the concept of transient networks. 



This theory was proposed in the seminal work [[EJ and developed in |4"6 , [4"7| , [48| . In 
the past decade, the theory of temporary networks was widely used to describe various 
observations from shear thickening of polymer solutions |^9[ to nonlinear relaxation of 



telechelic polymers ||50|| . The original theory was refined to account for (i) the effect of 
mechanical factors on the rates of breakage and reformation BTI, B31 KB] and (ii) the fi- 



nite rate of stress relaxation in dangling chains j54|. All previous studies, however, treated 
transient networks as consisting of chains with an identical number of strands. This hy- 
pothesis substantially simplifies the analysis, but it imposes unreasonable limitations on 
the model. Molecular dynamics simulation |55j reveals that the chain length dramatically 
affects its relaxation rate. Based on this observations, we assume a transient network to 
be composed of long chains with different lengths and characterize the network by the 
probability density of chains with various numbers of strands. This refinement provides 
a way to distinguish changes in the viscoelastic response driven by (i) an increase in the 
elastic modulus of a strand and (ii) the growth of the number of temporary crosslinks. 

The exposition is organized as follows. Section 2 deals with kinetic equations for the 
rates of breakage and reformation. The strain energy density of a transient network is 
determined in Section 3. Section 4 is concerned with constitutive equations for a tempo- 
rary network of long chains. The stress-strain relations are applied to describe uniaxial 
extension of a specimen in Section 5. Several hypotheses regarding (i) the distribution 
of chains with various number of strands and (ii) the effect of the number of strands on 
the rate of a chain's breakage are introduced in Section 6. The stress-strain relations are 
verified in Section 7 by comparing results of numerical simulation with experimental data 
for carbon black filled rubber in static and dynamic tests. Some concluding remarks are 
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formulated in Section 8. 



2 The concept of temporary networks 

A rubbery polymer is modeled as a network of long chains (with various numbers of 
strands) bridged by temporary junctions. A chain whose ends are connected to separate 
junctions is treated as an active one. Snapping of an end of a chain from a junction 
is thought of as its breakage (transition from its active state to the dangling state). 
When a dangling chain captures a junction, a new active chain arises [4~8"|. Breakage and 
reformation of active chains are treated as thermally activated processes: attachment and 
detachment events occur at random times as they are driven by thermal fluctuations. 

A chain is determined by two parameters: the number of strands, n = 1,2..., and 
the end-to-end vector, R. The distribution of chains with various numbers of strands 
is determined by the probability density p(n), which is assumed to be independent of 
temperature, T, and the strain intensity. This hypothesis is equivalent to the assumption 
that chains do not fall to pieces under loading. 

The kinetics of reformation is described by the function X(t,T,n), which equals the 
number of active (merged to the network) chains with n strands (per unit mass) at time 
t that has last been linked to the network at instant r G [0, £]. This function entirely 
determines the current state of a rubbery polymer at the micro-level. For example, the 
quantity X(t,t,n) equals the number of active chains (per unit mass) with n strands at 
time t. In particular, X(0, 0, n) is the initial concentration of active chains with n strands, 

X(0,0,n) = Ep(n), (1) 

where H is the total number of active chains per unit mass (we suppose that S does 
not change with strains and coincides with the number of active chains in a stress-free 
medium). The amount 

— (t,T,n) dr 
or t =T 

is the number of active chains with n strands (per unit mass) that have been linked to 
the network within the interval [r, r + dr] ; the quantity 

dX , 

— (t,r,n) dr 

OT 

determines the number of these chains that did not break during the interval [r, t] ; the 
amount 

dX , 
- — (t,0,n) dt 

is the number of active chains (per unit mass) that detach from the network (for the first 
time) within the interval [t,t + dt], and the quantity 

9 2 X , 

-{t, t, n) dtdr 



dtdr 
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is the number of long chains (per unit mass) that has last been linked to the network 
within the interval [r, r + dr\ and leave the network (for the first time after merging) 
during the interval [t, t + dt}. 

The kinetics of evolution for a network is determined by the relative rate of breakage 
of active chains, T(t,n), and by the rate of merging of dangling links with the network, 
7(t,n). The assumption that T and 7 depend on the number of strands n distinguishes 
the present model from previous ones where all chains are treated as composed of an 
identical number of strands. The dependence of the rates of breakage and reformation on 
time reflects the effect of the current strain intensity on these parameters at time-varying 
loading. 

The relative rate of breakage, T, is defined as the ratio of the number of active chains 
broken per unit time to the total number of active chains, 



r(*,n) = - Y ,.\ ~ (t,0,n) 
X(t,0,n) ot 



F(t,n) 



dX. 



-1 P)2 



d 2 X 
dtdr 



(t,r,n). (2) 



In the general case, T depends on the instant r when a chain has been last bridged to 
the network (before its detachment) and on the guiding vector I (the unit vector whose 
direction coincides with the end-to-end vector of the chain at the time r). We neglect these 
dependencies, which is equivalent to the assumption that the network is homogeneous (the 
rate of breakage for a chain is independent of the instant of its connection to the network) 
and isotropic (the rate of breakage for a chain is independent of its direction at the instant 
of merging) . 



In agreement with conventional theories of temporary networks ||48|| , the rate of merg- 
ing, 7, is defined as the number of dangling chains (per unit mass) bridged to the network 
per unit time 

dX 

7(*»") = -^r(^ T > n ) _ ■ ( 3 ) 

The hypothesis about the isotropicity of the network implies that 7 is independent of the 
guiding vector I. 

We consider a model with two states (active and dangling) of a long chain. A dangling 
chain is treated as a chain where stresses totally relax after its detachment (no memory 
about previous deformations), whereas the active chain is a chain that preserves entire 
memory about the macro-strain at the instant of its connection with the network. One 
can hypothesize about an intermediate state between these two extremities, when the 
interval between the detachment event and the subsequent attachment to the network is 
too small for total relaxation of stresses, but is sufficiently large to ensure that a part of 



these stresses does relax [54|. This is equivalent to the assumption that a partial memory 
preserves in a chain about the history of its deformation. For the sake of simplicity, we 
exclude this case from the consideration. However, because preserving a partial memory 
about the loading history seems quite natural, we distinguish the rate of transition from 
the active state to the dangling state, T, and the attempt rate, T*, which is determined 
as the relative rate of slippage from sticky junctions for ends of active chains. It is 
assumed that substantially exceeds the rate of breakage, T, which means that majority 
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of active chains detach from the network and, immediately afterwards (before stresses 
totally relax), merge with some junctions. 

Equations (2) may be treated as ordinary differential equations for the function X. 
Integration of Eq. (2) with initial conditions (1) and (3) implies that 



X(t,0,n) 
dX 



dr 



(t,r, n) 



Ep(n) exp 
7(Y, n) exp 



— / r(s, n)ds 
Jo 

— J T(s, n)ds . 



(4) 



To exclude the rate of reformation 7 from Eq. (4), we postulate that the concentrations 
of active and dangling chains remain constant and independent of the strain level. This 
means that the number of chains broken per unit time coincides with the number of chains 
merging with the network within the same interval. The number of long chains detached 
from the network during the interval [t, t + dt] equals the sum of the number of initial 
chains (i.e., not broken until time t) that slip from temporary junctions 



dX 
~dt 



(t,0,n) dt 



and the number of chains connected with the network within the interval [r, r + dr] (for 
the last time before instant t) and broken within the interval [t, t + dt] 



d 2 X 



(t, t, n) dtdr. 



dtdr 

This number coincides with the number of chains attached to the network within the 
interval [t,t + dt], 

j(t, n) dt, 

which results in the balance law 



/x 9X , n , r l crx . 
7(t,n) = - — (t,0,n)-y o ^(t,r,n)dr. 



t d 2 X 



dt Jo dtdr 

Substitution of Eq. (4) into this equality implies that 



l{t,n) = T(t,n)\Ep(n)exp 



J T(s,n)ds + J 7(7-, n) exp — J T(s,n)ds drj. 



To solve this equation, we introduce the notation 



7(t, n) = 7(t, n) exp 



T(s, n)ds 



and find that 



Ep(n) + / 7(r, n)di 
Jo 



7(t,n) = T(t,n) 
Setting t = in Eq. (6), we obtain 

7(0,n) = r(0, n)Ep(n). 



(5) 
(6) 

(7) 
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Differentiation of Eq. (6) with respect to time results in 



, . 1 9T. . 
T(t,n) + — — — — (t,n) 



l(t,n). 



T(t,n) dt 

Integration of this equality from zero to t yields 

7(0, n) r(0,n) Jo 
It follows from Eqs. (7) and (8) that 



(8) 



/ r(s,n)ds 

j0 



7(£, n) = Sp(n)T(t, n) exp 

Combining this equality with Eq. (5), we arrive at the formula for the rate of reformation 

j(t,n)=Ep(n)F(t,n). (9) 

3 Strain energy density of a temporary network 

An active chain is thought of as a nonlinear elastic spring whose position is determined 
by the unit vector I directed along the end-to-end vector R. We suppose that the stress in 
a dangling chain totally relaxes before this chain captures a new junction and the natural 
(stress-free) configuration of a chain merging with the network at time r coincides with 
the actual configuration of the network at that instant. 

We begin with the calculation of the extension ratio (along the chain), A, for a chain 
that merged with the network at time r and has not been broken within the interval [r, t] . 
At instant r, the chain has a small length 5 and connects two junctions at points A and 
B. At time t > r, the junctions occupy points with the radius vectors 

fA(t) = f A {r) + u(t, r, f A (r)), f B (t) = f B (r) + u(t, r, f B (r)), 

where u(t, r, f) is the displacement vector at point f for transition of the network from the 
deformed configuration at time r to the deformed configuration at time t. The position 
of the chain is determined by the end-to-end vectors 



R(t) = f B (r)-f A (r)=5l, 

R(t) = f B (t)-f A (t)=5l+[u(t,T,f B (T))~u(t,r,f A (r))] 
= 61+ [u(t,r,r A (r) + 81) - u(t,r,r A (r))}. 



(10) 



Neglecting terms of the second order of smallness compared to 5, we find from Eq. (10) 
that 

R(t)=6l-[i + V(T)u(t,T,f A (T))\, 

where V(t) is the gradient operator in the actual configuration at instant t, I is the unit 
tensor and dot stands for the inner product. In terms of the radius vector f, this equality 
reads 



R(t) = Sl-V(r)f(t) = V(r)f(t) • 51, 
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where T stands for transpose and the argument ta is omitted for simplicity. The end-to- 
end length of the chain, ds, is given by 

ds 2 = R- R. 

Substitution of Eqs. (10) and (11) into this equality results in 

ds 2 {r) = 5 2 , ds 2 {t) = 5 2 l- C{t, t) ■ I, 



where 



C(t, 



V(r)f(t)- V(r)f(t) 



is the relative Cauchy deformation tensor for transition from the deformed configuration 
of the network at time r to the deformed configuration at time t. The extension ratio A 
is determined as 

Kt,r,T)- ds{t) 



ds(r) 



l-C{t,r)-l 



(12) 



With reference to IBS 



we suppose that the network is incompressible and neglect the 
energy of interaction between long chains. These hypotheses are based on the assumption 
that the effect of excluded volume and other multi-chain effects for an individual chain 
are screened by surrounded macromolecules. 

The mechanical energy of a network is defined as the sum of the mechanical energies 
for elongation of long chains comprising the network. This assumption is equivalent to the 
hypothesis that the bending energy of semi-flexible chains is small compared to its energy 
of elongation and may be neglected in the formula for the total energy of a transient 



network. It is worth noting some recent studies p7[ |58| where the bending energy is 
included into the formula for the mechanical energy of an ensemble of chains. 

The mechanical energy of a chain with n strands that merged with the network at 
instant r and has not been broken until the time t reads 



w(X{t,T,l),n 



where W is a smooth function. We postulate that W(X,n) may be factorized as 

W(X,n)=fi(n)W (X), (13) 

where /x(n) is the rigidity of a chain with n strands (with the dimension of energy) and 
Wo (A) is a dimensionless function of the extension ratio which satisfies the condition 



W Q (1) = 0. 



(14) 



Equation (14) means that the mechanical energy vanishes for a non-deformed chain. 

The mechanical energy (per unit mass) of initial chains which have not been broken 
within the interval [0,t] is given by 



1 oo „ 

— J2Kn)X(t,0,n) / W (X(t,0,l))dA(l), 

47F n=l JS 



(15) 
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where S is the unit sphere in the space of vectors / and dA{l) is the surface element on 
S. The mechanical energy (per unit mass) of long chains that merged with the network 
within the interval [r, r + dr] and have been connected with the network until the current 
time t reads 

i oo r)Y" r _ 

4^EM«)-^r(*,T,n)*- J s W (X(t,T,l))dA(l). (16) 

Summing expressions (15) and (16), we arrive at the formula for the strain energy density 
(per unit mass) of the network of long chains 



1 OO r „ 

U (t) = T-EM™) X(t,0,n) / W (X(t,0,l))dA(l) 
+ J o —(t,r,n)dT J s W (\(t,r,l))dA(l) . 
It follows from Eq. (12) that for an isotropic network, the expression 

W (X(t,r,l))dA(l) 



(17) 



depends on the principal invariants of the relative Cauchy deformation tensor, C(t,r), 
only. Because the third principal invariant of this tensor equals unity for an incompressible 
medium, we obtain 

/ W (X(t,T,l))dA(l) = 4irW(l k (t,r)), (18) 

where Ik(t,r) is the kth principal invariant of C(t,r) (k = 1,2), and W = J2) is 

the average (over end-to-end vectors) mechanical energy (at time t) of chains that merge 
with the network at instant r. Equations (13) and (18) imply that 

W(3,3) = 0, (19) 

where the conditions I± — 3 and I 2 — 3 describe the stress-free state of a chain. 
By analogy with Eq. (18), we write 

/ W (X(t,0,l))dA(l) =4nW(l 0k (tj), (20) 

where Iokit) is the kth principal invariant of the Cauchy deformation tensor for transition 
from the initial (stress-free) configuration of the network to its deformed configuration at 
time t. 

Substitution of Eqs. (18) and (20) into Eq. (17) results in the formula 



ft f)V 

U(t) = ]T »(n)\x(t,0,n)w(l ok (tj) + J q —(t,r,n)W (l k (t,r))dr 



(21) 



Our objective now is to calculate the derivative of the function U with respect to time. 
It follows from Eqs. (2), (19) and (21) that 

^(t) = Mt) - J 2 (t), (22) 



where 



and 



Ji(t) = £//(n)[x(t,0,n) — (I ofc (t)) + j Q —(t,T,n) — (I k (t,T))dr 



00 r pt f)X 

Mt) =^(n)r(t,n)\x(t,0,n)W(l ok (t)) + j Q —(t,r,n)W{l k (t,r))dr 



(23) 



(24) 



Bearing in mind that the principal invariants of the relative Cauchy deformation tensor 
C(t, r) coincide with the principal invariants of the Finger tensor for transition from the 
deformed configuration at time r to the deformed configuration at time t, 

F(t,r) = [V(r)f(t)] T - V(r)f(t), 
and using the chain rule for differentiation, we find that 



lr W«.r))=Efl^(4(«,r))^(F( t ,r)): 



OF, i J 



(25) 



where the colon stands for convolution of tensors. It is easy to prove that 

f(F)=/, f(F) = / l( F)/-^. 

Substituting these expressions into Eq. (25) and taking into account the symmetry of the 
Finger tensor, we obtain 



dW 



f dW r dW x t dW r 



OF 
~dt' 



dh " di 2 > di 2 

The derivative of the Finger tensor with respect to time is given by [[59 
dF 



(26) 



Of 



(V 



where 



V(t)v(t) .F(t,r)+F(t,T)-V(t)v(t), 



*(*) = §(*) 



(27) 



is the velocity vector. It follows from Eqs. (26) and (27) that 



dW 



dW _ dW 



fUVV 



dh 



dW 



: D, 



(2f 



where 



D(t) = - 
K ' 2 



V(t)v(t)) +V{t)v(t) 
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is the rate-of-strain tensor. Substitution of Eq. (28) into Eqs. (22) and (23) results in 



dU . s 



n=l 



J 2 (t) + 2 £ fi(n) X(t, 0, n) U i(*)F (t) + ^ 02 (t)F 2 (t) 



+ 



* dX 
o dr 



where 



— (t, r, n) [V>i(t, r)F(t, r) + ^ 2 (t, r)F 2 (i, r) 



«9iy r dw 



dr\:D(t), (29) 



(30) 



and the subscript "zero" stands for the quantities that describe transition from the stress- 
free configuration of the network (before loading) to its deformed configuration at the 
current time t. 



4 Constitutive equations for a transient network 

The absolute temperature T is assumed to be close to its reference value T , which implies 
that the effect of temperature on material parameters, as well as thermal expansion of 
the network may be neglected. For an incompressible network, the first law of thermody- 
namics reads [B0| 

where p is the mass density, <7<j is the deviatoric component of the Cauchy stress tensor a, 
q is the heat flux vector, $ is the internal energy and r is the heat supply per unit mass. 



The Clausius-Duhem inequality implies that [60 

dQ dS - /g\ />r 
p -dt =P lt +V <f)-T -°' (32) 

where S is the entropy and Q is the entropy production per unit mass. Bearing in mind 
that 

V.(|) = iv.g-ig.vT, 
and excluding the term V • q from Eqs. (31) and (32), we find that 

dQ dS d$ 1/ A ~ 1 -T7r r \^n /qq\ 

The internal energy $ is given by the conventional formula 

$ = ^ + ST, 

where ^ is the free (Helmholtz) energy per unit mass. Substitution of this equality into 
Eq. (33) yields 

dQ „dT dV 1/ A 1 - n , , 
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The following expression is accepted for the function 



q = V + (c-S )(T-T 



cTln^ + U, 



(35) 



where c is a specific heat, and an d S are the free energy and the entropy per unit 
mass in the stress- free configuration at the reference temperature T . Formula (35) means 
that the configurational entropy (associated with stretching of long chains) is neglected 
compared to mechanical energy and only the entropy induced by thermal motion 



S 



T 

S + cm — 

-10 



(36) 



is taken into account. The use of this hypothesis may be explained by the difficulties in 
calculating the increment of the configurational entropy driven by stretching of long chains 
§yj, as well as by some ambivalence in the definition of the mechanical energy, U, for 
nonlinear entropic springs . The latter means that one cannot distinguish contributions 
of the configurational entropy for non-Gaussian chains and their mechanical energy into 
Eq. (35) for the free energy per unit mass, unless physical assumptions are specified for 
the choice of the function W in Eq. (21). 

Substitution of expressions (29), (35) and (36) into Eq. (34) results in 



T dQ 
dt 



Mt) - 



P T(t) 



q(t) • V(t)T(t) 



+ -|<x d - 2p £ fi(n) X(t, 0, n) (i/) 01 (t)F (t) + il> m {t)F%(t) 



n=l 



t dX 



~ L ^(t,r,n)(^i(t,r)F(t,r)+^ 2 (t,r)F 2 (t,r))rfr 



: D(t). (37) 



It follows from Eqs. (24) and (37) that for an arbitrary loading program, the rate of 
entropy production is nonnegative, provided that 

1. the heat flux vector q obeys the Fourier law 

q(t) = -«V(t)T(t) 
with a positive thermal diffusivity k, 

2. the Cauchy stress tensor a satisfies the constitutive equation 



a(t) = -P(t)I + 2p£>(n) X(t,Q,n)(4> 01 (t)F (t) + i;o2(t)F2(t) 

rt 8X 



n=l 



+ 



dr 



(t, r, n) Mt, r)F(t, r) + ^ 2 (t, r)F 2 (t, r))dr 



(3* 



where P(t) is pressure. 
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To employ Eq. (38) as a constitutive equation for a rubbery polymer, one should establish 
connections between the relative Finger tensor for a temporary network, F(t, r), and that 
for the polymer at the macro-level. The simplest assumption regarding these tensors is the 
affinity hypothesis which postulates that F coincides with the Finger tensor for macro- 
strains in a specimen. A conventional explanation for this statement is that "surrounding 
molecules suppress the movements of the crosslinks so strongly that their positions change 
affinely with the shape of the specimen" . 

For a permanent network, where the rate of breakage and reformation of chains vanish, 
Eqs. (4) and (38) are transformed into the Finger formula for the Cauchy stress tensor 



a = -pi + 2aU 01 f + ^j 02 f^ 



with 



n 



(39) 



71=1 



5 Uniaxial tension of a specimen 

In this Section, Eq. (38) is employed to determine stresses in a sample at uniaxial exten- 
sion. Points of a polymeric specimen refer to a Cartesian frame {Xi} in the stress-free 
state and to a Cartesian frame {x{\ in the deformed state, (i = 1,2,3). Uniaxial tension 
of an incompressible medium is described by the formulas 



X\ = k(t)Xi, 



%2 



X2, %3 — f 7777 ) ^3, 



where k is the extension ratio. The relative deformation gradient V(r)f(i) reads 



^ (r)f(t) = wr 1 + im) te " 2e " 2+eV " 3 



where are base vectors of the frame {Xj}. The relative Finger tensor F(t,r) is given 
by 

- . . fk(t)\ 2 k(r) ( \ 

F(t,r)=(^) ei e 1 + ^(e 2 e 2 + e3e3). 

Substituting Eq. (40) into Eq. (38), we find the non-zero components 



(40) 



00 p 

P(i) +2p$>(n) X(t,0,n)(ij 01 (t)k 2 (t) +i/> 02 (t)k i 



n=l 



+ 



dr 



a 2 (t) 



dr v ' y V riV ' J \k(r)J rzv ' ' \fc(r) / 
00 r 

P(t) + 2p ]T n(n) X(t, 0, n) (^ m {t)k-\t) + ^{t)k~\t) 



(41) 
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of the Cauchy stress tensor 

a = (Tieiei + <x 2 (e 2 e 2 + e 3 e 3 ). 
Excluding the unknown pressure P(t) from Eq. (41) and the condition 

cr 2 (t) = 0, 

we arrive at the formula for the longitudinal stress 



(t) = 2pY / Kn)\x(t,0,n)[ij i(t)(k\t)-k-\t))+4> 02 (t)(k\t)-k- 2 (t) 



71=1 



+ 

+■ 



r*dX, , r , , ,ffk(t)\ 2 k(r)\ 



«-)((n) 4 -(lf) 



fc(t) 



(42) 



where the function X is determined by Eqs. (4) and (9) and the functions ■?/>/« are given 
by Eq. (30). 

We confine ourselves to the Mooney-Rivlin material with the mechanical energy 

m 

J 2 ) = Cl (/! - 3) + c 2 (J 2 - 3), (43) 

where ci and c 2 are adjustable parameters. The choice of the Mooney-Rivlin equation 
may be explained by two reasons: (i) formula (43) fits observations for some rubbery 
polymers and polymeric melts with a high level of accuracy |J5] and (ii) Eq. (43) provides 
a compromise between the concepts of affine and phantom networks |)6| . 
Substitution of Eqs. (30) and (43) into Eq. (42) results in 



a l {t) = 2pY,^{n) X{t^ 1 n)(c 1 + c 2 k~\t))(k 2 {t)-k-\t) 

71=1 

f l dX. . / 
+ / (t,T,n)(ci + c 2 



dr 



lo dr v "' ''-'V^ ' ~*k(t);\\k(T); k(t) 
Combining this equality with Eqs. (4) and (9), we arrive at the stress-strain relation 

oi(t) = 2pS^/x(n)p(n) exp(- / r(s, n)c?sj fcj + c 2 fc _1 (t)) (^ 2 (t) — A; -1 

n=l L JO 

-r r (-)-(-/>(-H(^^)((i) 2 -n 

For a permanent network, Eq. (44) is reduced to the conventional formula 

a x = 2(d + 02k- 1 ) {k 2 - kr 1 ) (45) 

with 



dr 



.(44) 



a = AC,;. 



(46) 
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We now analyze the viscoelastic response of a rubbery polymer, when small oscillations 
with a fixed amplitude k-y and a fixed frequency uj are superposed on the time-independent 
stretching of a specimen with the longitudinal elongation k , 



k(t) = ko + k\ exp(iu;t). 



(47) 



Assuming k to be of the order of unity and k x to be small compared to unity and 
neglecting terms of the second order of smallness, we find that 



ci + c 2 



k[r) 



kit) 
c l + c 2 k- l {t) 
( k(tU2 k(r) 



Cl + c 2 



1 — -^-(exp(Ljt) — exp(kjr)) 



ci + c 2 k 1 - ^c 2 exp(kut), 

Kq 

3^-(exp(k<jt) — exp(ioJT)) , 



~k(r); k(t) 
k 2 (t)-k^(t) 

It follows from these equalities that with the required level of accuracy, 
(ci + c 2 k-\t)) (k 2 {t) - k-\t)) = (ci + c 2 k^){k 2 - k, 1 ) 



{kl - k r ) + ^(2k 2 + k^ 1 ) exp(i^t). 
ko 



4 

k 



Cl {2k 2 + k ') + c 2 k Q - 1 {kl + 2k^) 



exp(io;t), 



/ k(r)\ / / k(t) \2 k(r)\ k\ , ,/ ,. . ,. 

Substitution of Eq. (48) into Eq. (44) implies that 

a 1 (t) = S (t) + S 1 (t) + S 2 (t), 

where 

oo 

S (t) = 2p( Cl + c 2 k r )(k 2 - k, 1 ) ]T t*(n)X(t, 0, n) 



(48) 



71=1 



determines relaxation of stresses in a specimen stretched with the extension ratio ko, 



Si(t) = 2p 



h 
ko 



( Cl {2k 2 + k^ 1 ) + c 2 ko\k 2 + 2&Q- 1 )) - 3 



x exp(ic;t) ^ p(n)X(t, 0, n) 



n=l 



characterizes relaxation of small oscillatory stresses and 



u 00 
S 2 (t) = 6p-i(ci + c 2 )^/iH 

fc n= i 

r* dX 



X(t, t, n) exp(iu;i) 



P OX . , .. , , 

— / — — It, t, n) exp{iujT)dT 
Jo or 

ki 00 
6 Py-( c i + c 2)2 J! ^(n)p(n) 
fc o „=i 

— y T(r, n)exp^— y T(s, n)rfs — icj(t — r)^c?r 



exp(ia;t) 
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determines the steady response to small oscillations superposed on uniaxial stretching. 
Following f23fl , we define the complex elastic modulus E* as the ratio of S 2 (t) to the 
oscillatory strain 

= T-exp(ic<jt) 

and obtain 



1 



E*(t,uj) = 6p(ci + c 2 )EJ2 K n )p( n 

n=l 

— J r(r, n)exp^— J r(s,n)ds — iu(t — r)^dr 



(49) 



For the loading process (47), the rate of breakage V is time- independent, T = T(n), (it is 
determined by the time-independent "basic" elongation ko)- Introducing the new variable 
s = t — r and replacing the upper limit of integration in Eq. (49) by infinity (in agreement 
with the conventional procedure for the analysis of dynamic response), we obtain 

r c oo 

E*(uj) = 6p( Cl + c 2 )EJ2 K n )p( n ) l-T(n) / exp(-(r(n)+iw)(t-r))dr 

n=l Ju " ' 

oo 

= 6p(ci + c 2 )S K n )p{ n 



n=l 



1 - 



Tin) 



Tin) 



1LU 



We split the complex modulus E* into the sum of storage and loss moduli, 

E*(lu) = E'(uj)+\E"(uj), 

and arrive at the formulas 

^ fi(n)p(n)uj 2 



E\u) 
E"(u) 



6p(ci + c 2 )H J2 

n=l 

oo 

6p(ci + c 2 )S 



fi(n)p(n)T[n)LJ 



71=1 



r 2 (n) +u 2 



(50) 



It follows from Eq. (50) that for a Mooney-Rivlin viscoelastic medium, the dynamic 
moduli are independent of the elongation, ko, on which small oscillations are superposed. 
This result is a substantial shortcoming of Eq. (43) because experimental data reveal that 
the storage and loss moduli noticeable change with stretching. 



6 Adjustable parameters in the model 

Governing equations (45) and (50) are determined by two adjustable parameters, C\ and c 2 , 
which characterize the elastic response of a polymer at finite strains, and three material 
functions, pin), fi{n) and T{t,n). The function p{n) describes the distribution of long 
chains with various numbers of strands. This function is analogous to the probability 
density of traps with various energies in the energy-landscape theory for structural glasses 
Because the distribution of long chains in rubbery polymers is determined by the 
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vulcanization process (whose description at the micro-level is far from being exhausted, 
see |67| and the references therein), no explicit expression has been yet proposed for this 



function. As a first approximation, we assume the distribution of chains to be exponential 

exp(-an) 

P{n) = t-^ r, (51) 

exp(a) — 1 

where a is an positive constant. The factor exp(a) — 1 in Eq. (51) is determined from 
the condition 

oo 

X>0) = i 

n=l 

and the formula 

00 00 n 1 1 

exp(-an) = ^ [exp(-a)] n - 1 = - 1 = — -. (52) 

^o L J l-exp(-a) exp(a)-l 

The average number of strands in a chain is given by 

oo j oo 1 (9 °° 

( n ) = np(n) = — J2 n exp(-an) = — -— ^ exp(-an). 

exp(a) - 1 ^ exp(a) - 1 da ^ 

Combining this equality with Eq. (52), we find that 



< 53 > 

The other adjustable function, fi(n), characterizes the longitudinal rigidity of a long 
chain containing n strands. With reference to the Rouse model, a chain may be treated 
as a system of n identical springs connected in sequence. The rigidity of such a system is 
given by the conventional formula 

Ai(n) = ^, (54) 
n 

where /xq is the rigidity of an individual spring. We suppose that Eq. (54) may be applied 
(as a first approximation) to long chains with an arbitrary geometry. 

The third function to be found is the rate of relaxation F(t, n). It is assumed that 
r(t,n) is factorized as 

T(t,n) = T (t)r ] (n), (55) 

where To is responsible for mechanically-induced changes in the relaxation rate, whereas 
77(71) describes the effect of the chain's length on the rate of its slippage from a temporary 
junction. We suppose that rj exponentially grows with the number of strands n, 

77(71) = exp(/?n), (56) 

where (3 is a positive parameter. Equation (56) is in contrast with the Eyring formula 
which is conventionally adopted for the description of thermally activated processes in 
polymers [68]. Within the concept of traps, the rate of relaxation is traditionally presumed 
to decrease exponentially with the energy of a potential well on the energy landscape where 
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a relaxing region is trapped. Assuming this energy to be proportional to the volume of 
a micro-domain, we find that the relaxation rate decreases with the number of strands 
participating in collective rearrangement, whereas Eq. (56) implies that this rate grows 
with the number of strands in a long chain. This contradiction may be explained with 



reference to the reptation theory for polymeric chains [p9 |. A semiflexible chain with a 
finite bending stiffness is thought of as a curvilinear rod whose micro-motion is restricted 
to a tube composed by surrounding chains. The radius of the tube is estimated as a few 
lengths of a strand, which implies that the micro-motion of a chain may be split into 
diffusion along the tube and small lateral fluctuations around the tube's centerline. For 
a stiffless rod, lateral fluctuations do not affect the response of junctions, which implies 
that breakage of chains is mainly associated with longitudinal diffusion. In this case, 
small thermal fluctuations induce rather small longitudinal displacements whose energy 
is insufficient to detach a chain with a large number of strands from the junctions. The 
opposite picture is observed for a chain with a finite bending stiffness, whose lateral 
oscillations (driven by relatively small, but random fluctuations) may be amplified due to 
interaction of transverse elastic waves in the rod and may cause tearing of a chain from 
temporary junctions. This study deals with long chains with a non-vanishing bending 
stiffness, for which Eq. (56) serves as a natural hypothesis. 

Substitution of expressions (51) and (54) into Eq. (39) results in 



Bearing in mind that 



_ p p5 ~ exp(-cm) 
exp(a) - 1 n 

exp(—nx)dx, 



exp(— an 



oo 



n 

we rewrite this equality as follows: 

A _ PopS ( p pH [°° dx 



expect J 



/ > exp(— ax)dx = — — / — . 

— 1 J a ~y exp(a) - 1 ia exp(:r) — 1 ' 



where the sum is determined with the help of Eq. (52). To calculate the integral, we 
introduce the new variable y = exp(x) and obtain 

f°° dx f°° dy f 00 ( 1 K , , / . . \ 

/ T~\ = / ~i 7 = / ( lay = a — ln(exp(aj — 1). 

Ja exp(x) — 1 Jexp(a) y(y — 1) Jexp(a) K y — 1 y' v ' 

This implies that 

„q-ln(exp(a) -1) 

A = poP^ -r-, : • 57 

exp(a) — 1 

Setting 

" — 

where S is the number of strands per unit mass, we find that 

A = Ao -ln(exp(a)-l) 
(n)(exp(a) — 1) 
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with A = /UopS . Equations (53) and (58) allow the average rigidity of an ensemble 
of chains, A, to be expressed as a function of the average number of strands in a chain 
(n). Figure 1 reveals that the quantity A [an analog of the elastic modulus of a rubbery 
polymer, see Eq. (46)] decreases with the average number of strands in a chain (n). With 
an acceptable level of accuracy this dependence may be described by the scaling law 

A oc (n)- K 

with k = 0.56. 

7 Validation of the model 

Our objective now is to find adjustable parameters of the model by fitting observations for 
a carbon black filled rubber at various temperatures, T, and various programs of loading. 

7.1 Tensile tests with high rates of loading 

We begin with matching observations in tensile tests with relatively high rates of loading 
(when the processes of breakage and reformation of long chains during the test may be 
neglected) and moderate deformation (when mechanically induced alignment of chains is 
not observed and the network may be treated as isotropic). For a detailed description of 
specimens and the experimental procedure, see [ |T6|j . 

The loading process is characterized by the rate of engineering strain 

. _ de eng _ 

c eng > c eng ^ ± - 

Experimental data in tensile tests with the strain rate e eng = 0.2 s _1 at four different 
temperatures are plotted in Figure 2 together with the results of numerical simulation. 
The parameters C\ and C2 in Eq. (45) are determined using the least-squares algorithm. 

It is found that the quantity C\ may be set to be zero, and the model with only one 
adjustable parameter, C2, correctly predicts observations at temperatures from T = 296 K 
to T = 373 K. At the lowest temperature, T = 253 K, experimental data slightly deviate 
from the model prediction. These discrepancies may be explained by (partial) failure of 
carbon filled rubber and insufficient accuracy of the Mooney-Rivlin equation (43). 

The influence of temperature T on the elastic modulus is illustrated by Figure 3, where 
the quantity C2 is plotted versus the degree of undercooling AT = T — T g . The elastic 
modulus monotonically decreases with temperature until some critical temperature T cr 
and, afterwards, remains practically constant. In the interval [T g ,T cr ] the dependence 
Ct,(T) is fairly well approximated by the linear function 

C 2 = a - ax AT, (59) 

where the quantities and a\ are found using the least-squares technique. The decrease 
in C2 with temperature confirms our hypothesis that in the region of temperatures under 
consideration, the effect of configurational entropy (whose contribution into the free energy 
linearly increases with temperature) is negligible compared to the mechanical energy of 
chains. 
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7.2 Dynamical tests on a compressed specimen 

As it was discussed in Introduction, an increase in the elastic modulus of a polymer with 
a decrease in temperature is traditionally associated with two phenomena: (i) an increase 
in the rigidity of a strand, /j,q, and (ii) an increase in the number of temporary junctions 
(which is tantamount to a decrease in the average number of strands in a chain, (n)). 
To assess which mechanism is responsible for temperature-dependent changes in elastic 
moduli and to analyze the effect of temperature on the rate of breakage of chains, we 
fit experimental data for a preloaded specimen (the longitudinal strain eo = —0.1) in 
dynamic tests with the small amplitude of oscillations Ae = 0.006. A detailed description 
of the experimental procedure can be found in |[23|| . 

Substitution of Eqs. (51), (54) and (56) into Eq. (50) results in 

el// \ u 2 exp(-an) 

E (uj = C 2^ — — — — -, 60 

^ n(r exp(/3ra) + uj 2 ) 

where 

67r/! pH 

C = -— -ci + c 2 . 61 

exp(a) — 1 

Comparing Eqs. (46) and (61) and bearing in mind Eq. (57) and the equality C\ = 0, we 
find that 

a — ln(exp(a) — 1) 

Equation (62) provides a simple relation between the moduli observed in a uniaxial tensile 
test with a constant rate of strain and in a dynamic test with a small amplitude of 
oscillations. However, this formula should be taken with caution, because its derivation 
is based on the assumption about the Mooney-Rivlin equation for the strain energy 
density of a chain (the hypothesis which does not guarantee an acceptable quality of 
fitting observations in tensile tests with a constant rate of strain, and which may lead to 
even larger discrepancies between the linearized equation (50) and experimental data in 
dynamic tests). It is easy to show that for any strain energy density of a chain, W, [not 
necessarily determined by Eq. (43)], the parameter C is determined by Eq. (61), where 
the sum c\ + c 2 is replaced by some coefficient, which, in general, depends on the strain 
eo, but which is independent of temperature, T, and the average number of strands, (n). 

To check which of the two assumptions about the influence of temperature is more 
adequate for the description of the response of rubbery polymers, we suppose that the 
number of strands is temperature-independent (which is tantamount to the postulate 
that a is independent of temperature) and fit observations in dynamic tests at various 
temperatures. We set a = 0.02, which corresponds to the average number of strands 
(n) = 1.24 x 10 5 (the value which is in agreement with available experimental data for 
rubbery polymers flSSH). The parameter (3 = 2.27 is found by fitting data in a test at 
ambient temperature by using the steepest-descent procedure. We fix this value of j3 and 
approximate observations at other temperatures with only two adjustable constants, C 
and r . Given a parameter C, the relaxation rate T is determined using the steepest- 
descent algorithm. An analog of rigidity, C, is found by the least-squares technique. 
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Figure 4 demonstrates good agreement between experimental data and results of nu- 
merical simulation with a time-independent parameter a. This means that our assump- 
tion that the number of strands remains constant at all temperatures in the region under 
consideration is quite acceptable for matching observations. 

The parameter C is plotted versus temperature in Figure 3 which demonstrates that 
in the region of temperatures [T g ,T cr ], the dependence C(T) is fairly well approximated 
by the linear function 

C = b - feiAT, (63) 

where the parameters bo and b\ are determined by the least-squares technique. Because 
the slopes of the graphs C^AT) and C(AT) are close to each other (the ratio ai/ao is 
0.0071, whereas 61/60 equals 0.0067), we may conclude that in the interval [T g ,T cr ], the 
effect of temperature on a is rather weak. This means that the average number of strands 
in a chain, (n), feebly depends on temperature, which implies that entanglements are 
practically not transformed into temporary junctions. 

The relaxation rate r is plotted versus the degree of undercooling AT in Figure 5. 
Experimental data are fairly well approximated by the function 

\ogT = do-d 1 AT, (64) 

where adjustable parameters do and d\ are found using the least-squares algorithm. 

Equation (64) implies that the rate of breakage for long chains decreases with tem- 
perature (in contrast with the theory of thermally activated processes which predicts the 
growth of r ). This result may be explained by a decrease in the bending stiffness of 
chains with temperature. At low temperatures, when the bending rigidity is relatively 
high, any local thermal fluctuation produces oscillations of a chain (which is thought 
of as a curvilinear elastic rod). Interaction of transverse oscillations driven at random 
times at random points of the rod with fixed ends amplifies their amplitude, and, as 
a consequence, induces relatively large displacements of the chain's ends (local thermal 
fluctuations are transformed into global ones at the length-scale of a chain). The ampli- 
fication of random displacements of the chain's ends may be sufficient for their slippage 
from temporary junctions, which is reflected in the high rate of breakage, r , in the vicin- 
ity of the glass transition point. With the growth of temperature, the bending rigidity 
substantially decreases, which results in a decline of transverse oscillations of chains. Lo- 
cal thermal fluctuations weakly interact and their amplitudes are not amplified. As a 
result, only thermal fluctuations in the close vicinity of a chain's ends can induce their 
slippage from temporary junctions. Although the number of these fluctuations increases 
and the strength of junctions decreases with temperature, the total number of broken 
chains (per unit time) diminishes, in agreement with data depicted in Figure 5. 

8 Concluding remarks 

Constitutive equations are derived for the nonlinear viscoelastic response of rubbery poly- 
mers at finite strains. The model is based on the concept of temporary networks, which 
treats an amorphous polymer as an ensemble of long chains connected to junctions. At 
random times, active chains detach from the junctions as they are thermally agitated. A 



21 



dangling chain merges with the network when its free end captures some junction in its 
vicinity. 

Unlike conventional concepts of transient networks, we suppose that (i) long chains 
consist of various numbers of strands and (ii) the rigidity of a chain and the rate of its 
breakage substantially depend on the chain's length. Several hypotheses are introduced 
regarding the distribution of chains with various lengths and the dependence of the rate 
of reformation on the number of strands. These assertions are verified by comparison 
with observations in uniaxial tensile tests. Fair agreement is demonstrated between ex- 
perimental data for a carbon black filled rubber at various temperatures and results of 
numerical simulation. 

The following conclusions are drawn: 

1. The average number of strands in a chain weakly depends on temperature. This 
implies that changes in the viscoelastic response of rubber with temperature may 
be ascribed to the effect of temperature on the average rigidity of a strand only. An 
acceptable agreement with experimental data is achieved without the hypothesis 
about a substantial increase in the number of junction with a decrease in tempera- 
ture (induced by transition of entanglements into temporary junctions). 

2. Some critical temperature, T cr , is found above the glass transition point T g . In the 
region [T g ,T cr ], the rigidity of a strand linearly decreases with temperature, while 
above the critical point, this parameter is independent of temperature. The values 
of T cr determined in static and dynamic tests are close to one another. 

3. Conventional theories of rubber elasticity are based on the assumption about the 
entropic nature of free energy. This implies that stresses linearly increase with 
temperature. Because this assertion contradicts experimental data for carbon black 
filled rubber, it is postulated that the entropic contribution into the free energy 
is small compared to that for the mechanical energy. This hypothesis is fairly well 
confirmed by observations, which demonstrate a decrease in stresses with the growth 
of temperature. 

4. Unlike the theory of cooperative relaxation in glassy polymers which presumes that 
the rate of rearrangement exponentially decreases with volume of a rearranged do- 
main, the rate of breakage for rubbery polymers increases with the growth of the 
number of strands in a long chain. This may be explained by a finite bending 
rigidity of a chain which is not taken into account by conventional theories. At 
relatively large bending rigidity, an active chain may be thought of as a curvilinear 
elastic rod. Local thermal fluctuations result in transverse oscillations of the rod. 
The amplitude of oscillations grows because of their interaction, which results in an 
increase in transverse displacements. The growth of the bending rigidity and the 
rod's length leads to an increase in the amplitude of oscillations for the rod's ends, 
and, as a consequence, to an increase in the rate of breakage. 

5. In contrast with the theory of thermally activated processes, the rate of breakage for 
active chains decreases with temperature. This may be explained by a substantial 
decline in the bending rigidity of chains which implies that local random oscillations 
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do not interact, and only thermal fluctuations in the close vicinity of the end points 
cause detachment of an active chain from the network. 

Acknowledgement 

AD gratefully acknowledge financial support by the Israel Ministry of Science through 
grant 1202-1-00. 



23 



References 



Govinjee, S., Simo, J.: A micro-mechanically based continuum damage model for 
carbon black-filled rubbers incorporating Mullins effect. J. Mech. Phys. Solids 39, 
87-112 (1991). 

So, H., Chen, U.D.: A nonlinear mechanical model for solid-filled rubbers. Polym. 
Eng. Set. 31, 410-416 (1991). 

Govinjee, S., Simo, J.: Mullins' effect and the strain amplitude dependence of the 
storage modulus. Int. J. Solids Structures 29, 1737-1751 (1992). 

Arruda, E.M., Boyce, M.C.: A three dimensional constitutive model for the large 
stretch behavior of rubber elastic materials. J. Mech. Phys. Solids 41, 389-412 (1993). 

Johnson, M.A., Beatty, M.F.: The Mullins effect in uniaxial extension and its influ- 
ence on the transverse vibration of a rubber string. Continuum Mech. Thermodyn. 
5, 83-115 (1993). 

Johnson, M.A., Beatty, M.F.: A constitutive equation for the Mullins effect in stress 
controlled uniaxial extension experiments. Continuum Mech. Thermodyn. 5, 301-318 
(1993). 

Lee, B.J., Argon, A.S., Parks, D.M., Ahzi, S., Bartczak, Z.: Simulation of large strain 
plastic deformation and texture evolution in high density polyethylene. Polymer, 34, 
3555-3575 (1993). 

Mukhopadhyay, K., Tripathy, D.K., De, S.K.: Dynamic mechanical properties of 
silica-filled ethylene vinyl acetate rubber. J. Appl. Polym. Sci. 48, 1089-1103 (1993). 

Le Tallec, P., Rahier, C, Kaiss, A.: Three-dimensional incompressible viscoelastic- 
ity in large strains: formulation and numerical approximation. Comp. Meths. Appl. 
Mech. Engng. 109, 233-258 (1993). 

Witten, T.A., Rubinstein, M., Colby, R.H.: Reinforcement of rubber by fractal ag- 
gregates. J. Phys. E France 3, 367-383 (1993). 

van den Bogert, P.A.J. , de Borst, R.: On the behavior of rubberlike materials in 
compression and shear. Arch. Appl. Mech. 64, 136-146 (1994). 

Hausler, K., Sayir, M.B.: Nonlinear viscoelastic response of carbon black reinforced 
rubber derived from moderately large deformations in torsion. J. Mech. Phys. Solids 
43, 295-318 (1995). 

Ulmer, J.C.: Strain dependence of dynamic mechanical properties of carbon black- 
filled rubber compounds. Rubber Chem. Technol. 69, 15-47 (1995). 

Holzapfel, G., Simo, J.: A new viscoelastic constitutive model for continuous media 
at finite thermomechanical changes. Int. J. Solids Structures 33, 3019-3034 (1996). 



24 



[15] Huber, G., Vilgis, T.A., Heinrich, G.: Universal properties in the dynamic deforma- 
tion of filled rubbers. J. Phys.: Condens. Matter 8, L409-L412 (1996). 

[16] Lion, A.: A constitutive model for carbon black filled rubber: experimental investi- 
gation and mathematical representation. Continuum Mech. Thermodyn. 8, 153-169 
(1996). 

[17] Lion, A.: A physically based method to represent the thermo- mechanical behaviour 
of elastomers. Acta Mech. 123, 1-25 (1997). 

[18] Lion, A.: On the large deformation behaviour of reinforced rubber at different tem- 
peratures. J. Mech. Phys. Solids 45, 1805-1834 (1997). 

[19] Reese, S., Wriggers, P.: A material model for rubber-like polymers exhibiting plas- 
tic deformation: computational aspects and comparison with experimental results. 
Comp. Meths. Appl. Mech. Engng. 148, 279-298 (1997). 

[20] Spathis, G.: Non-linear constitutive equations for viscoelastic behaviour of elastomers 
at large deformations. Polym. Gels Networks 5, 55-68 (1997). 

[21] Zaroulis, J.S., Boyce, M.C.: Temperature, strain rate, and strain state dependence of 
the evolution in mechanical behaviour and structure of poly(ethylene terephthalate) 
with finite strain deformation. Polymer 38, 1303-1315 (1997). 

[22] Kaliske, M., Rothert, H.: Constitutive approach to rate-independent properties of 
filled elastomers. Int. J. Solids Structures 35, 2057-2071 (1998). 

[23] Lion, A.: Thixotropic behaviour of rubber under dynamic loading histories: experi- 
ments and theory. J. Mech. Phys. Solids 46, 895-930 (1998). 

[24] Reese, S., Govidjee, S.: Theoretical and numerical aspects in the thermo-viscoelastic 
material behaviour of rubber-like polymers. Mech. Time- Dependent Mater. 1, 357- 
396 (1998). 

[25] Septanika, E.G., Ernst, L.J.: Application of the network alteration theory for the 
modeling the time-dependent constitutive behavior of rubbers. 1. General theory. 
Mech. Mater. 30, 253-263 (1998). 

[26] Septanika, E.G., Ernst, L.J.: Application of the network alteration theory for the 
modeling the time- dependent constitutive behavior of rubbers. 2. Experimental ver- 
ification. Mech. Mater. 30, 255-273 (1998). 

[27] Liana, P.G., Boyce, M.C.: Finite strain behavior of poly (ethylene terephthalate) 
above the glass transition temperature. Polymer 40, 6729-6751 (1999). 

[28] Boyce, M.C., Socrate, S., Liana, P.G.: Constitutive model for the finite deforma- 
tion stress-strain behavior of poly(ethylene terephthalate) above the glass transition. 
Polymer 41, 2183-2201 (2000). 



25 



[29] Tzika, P.A., Boyce, M.C., Parks, D.M.: Micromechanics of deformation in particle- 
toughened polyamides. J. Mech. Phys. Solids 48, 1893-1929 (2000). 

[30] Kramer, H.H., Berger, L.L.: Fundamental processes of craze growth and fracture. 
Adv. Polym. Sci. 91-92, 1-68 (1990). 

[31] Steenbrink, A.C., van der Giessen, E.: On cavitation, post-cavitation and yield in 
amorphous polymer-rubber blends. J. Mech. Phys. Solids 47, 843-876 (1999). 

[32] Dorfmann, A., Burtscher, S.L.: Aspects of cavitation damage in seismic bearings. J. 
Struct. Engng. 126, 573-579 (2000). 

[33] Hobeika, S., Men, Y., Strobl, G.: Temperature and strain rate independence of 
critical strains in polyethylene and poly (ethylene- co- vinyl acetate). Macromolecules 
(in press). 

[34] Gaucher-Miri, V., Seguela, R.: Tensile yield of polyethylene and related copolymers: 
mechanical and structural evidences of two thermally activated processes. Macro- 
molecules 30, 1158-1167 (1997). 

[35] Hiss, R., Hobeika, S., Lynn, C, Strobl, G.: Network stretching, slip processes and 
fragmentation of crystallites during uniaxial drawing of polyethylene and related 
copolymers. A comparative study. Macromolecules 32, 4390-4403 (1998). 

[36] Wu, P.D., van der Giessen, E.: On improved network models for rubbery elasticity 
and their applications to orientational hardening in glassy polymers. J. Mech. Phys. 
Solids 41, 427-456 (1993). 

[37] Treloar, L.R.G.: The Physics of Rubber Elasticity. Oxford: Clarendon Press 1975. 

[38] Ferry, J.D.: Viscoelastic Properties of Polymers. New- York: Wiley 1980. 

[39] Doi, M., Edwards, S.F.: The Theory of Polymer Dynamics. Oxford: Oxford Univer- 
sity Press 1986. 

[40] Adam, G., Gibbs, J.H.: On the temperature dependence of cooperative relaxation 
properties in glass-forming liquids. J. Chem. Phys. 43, 139-146 (1965). 

[41] Dyre, J.C.: Energy master equation: a low temperature approach to Bassler's 
random-walk model. Phys. Rev. B 51, 12276-12294 (1995). 

[42] Monthus, C, Bouchaud, J. -P.: Models of traps and glass phenomenology. J. Phys. 
A: Math. Gen. 29, 3847-3869 (1996). 

[43] Sollich, P.: Rheological constitutive equation for a model of soft glassy materials. 
Phys. Rev. E 58, 738-759 (1998). 

[44] Struik, L.C.E.: Internal Stresses, Dimensional Instabilities and Molecular Orienta- 
tions in Plastics. Chichester: Wiley 1990. 



26 



[45] Green, M.S., Tobolsky, A.V.: A new approach to the theory of relaxing polymeric 
media. J. Chem. Phys. 14, 80-92 (1946). 

[46] Yamamoto, M.: The visco-elastic properties of network structure. 1. General formal- 
ism. J. Phys. Soc. Japan 11, 413-421 (1956). 

[47] Lodge, A.S.: Constitutive equations from molecular network theories for polymer 
solutions. Rheol. Acta 7, 379-392 (1968). 

[48] Tanaka, F., Edwards, S.F.: Viscoelastic properties of physically cross-linked net- 
works. Transient network theory. Macromolecules 25, 1516-1523 (1992). 

[49] Wang, S.-Q.: Transient network theory for shear-thickening fluids and physically 
cross-linked systems. Macromolecules 25, 7003-7010 (1992). 

[50] Serero, Y., Jacobsen, V., Berret, J.-F., May, R.: Evidence of nonlinear chain stretch- 
ing in the rheology of transient networks. Macromolecules 33, 1841-1847 (2000). 

[51] Phan-Thien, N., Tanner, R.I.: A new constitutive equation derived from network 
theory. J. N on- Newtonian Fluid Mech. 2, 353-365 (1977). 

[52] Fuller, G.G., Leal, L.G.: Network models of concentrated polymer solutions derived 
from the Yamamoto network theory. J. Polym. Sci.: Polym. Phys. Ed. 19, 531-555 
(1981). 

[53] Petruccione, F., Biller, P.: Rheological properties of network models with confi- 
guration-dependent creation and loss rates. Rheol. Acta 27, 557-560 (1988). 

[54] Palierne, J.-F.: Sticky dumbbells: from Hookean dumbbells to transient network. 
Rheol. Acta 36, 534-543 (1997). 

[55] Barsky, S., Slater, G.W.: A nonequilibrium molecular dynamic simulation of the 
time-dependent orientational coupling between long and short chains in a bimodal 
polymer melt upon uniaxial stretching. Macromolecules 32, 6348-6358 (1999). 

[56] Ahn, K.H., Osaki, K.: Mechanism of thear thickening investigated by a network 
model. J. Non- Newtonian Fluid Mech. 56, 267-288 (1995). 

[57] Morse, D.C.: Viscoelasticity of concentrated isotropic solutions of semiflexible poly- 
mers. 1. Model and stress tensor. Macromolecules 31, 7030-7043 (1998). 

[58] Morse, D.C.: Viscoelasticity of concentrated isotropic solutions of semiflexible poly- 
mers. 2. Linear response. Macromolecules 31, 7044-7067 (1998). 

[59] Drozdov, A.D.: Finite Elasticity and Viscoelasticity. Singapore: World Scientific 
1996. 

[60] Coleman, B.D., Gurtin, M.E.: Thermodynamics with internal state variables. J. 
Chem. Phys. 47, 597-613 (1967). 



27 



[61] Altenberger, A.R., Dahler, J.S.: Statistical mechanics of rubber elasticity. J. Chem. 
Phys. 92, 3100-3111 (1990). 

[62] Everaers, R.: Elasticity of c*-gels. J. Phys. E France 5, 1491-1500 (1995). 

[63] Mooney, M.: A theory of large elastic deformation. J. Appl. Phys. 11, 582-592 (1940). 

[64] Rivlin, R.S., Saunders, D.W.: Large elastic deformations of isotropic materials. 7. 
Experiments on the deformation of rubber. Phys. Trans. Roy. Soc. London 243, 
251-288 (1951). 

[65] Zang, Y.H., Muller, R., Froelich, D.: Interpretation of the rheological behaviour in 
elongation of uncrosslinked polystyrene melts in terms of the Mooney-Rivlin equa- 
tion. Polymer 27, 61-65 (1986). 

[66] Ngai, K.L., Roland, CM.: Junction dynamics and the elasticity of networks. Macro- 
molecules 27, 2454-2459 (1994). 

[67] Goldbart, P.M., Castillo, H.E., Zippelius, A.: Randomly crosslinked macromolecular 
systems: vulcanization transition to and properties of the amorphous solid state. 
Adv. Phys. 45, 393-468 (1996). 

[68] Krausz, A.S., Eyring, H.: Deformation Kinetics. New York: Wiley, 1975. 



28 




Figure 2: The longitudinal stress <j\ MPa versus the engineering strain e eng for carbon 
black filled rubber in tensile tests at various temperatures T K. Circles: experimental 



data [I8|. Solid lines: their approximation by Eq. (45). Curve 1: T — 253; curve 2: 



T = 296; curve 3: T = 333; curve 4: T = 373 
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Figure 3: The coefficients C2 MPa (unfilled circles) and C MPa (filled circles) versus 
the degree of undercooling AT K for a carbon black filled rubber. Circles: treatment of 
observations [|18|, ^3f. Solid lines: approximation of the experimental data by Eqs. (59) 
and 63). Curve la: a = 1.8097, a x = 0.0128; curve 2a: b = 462.16, 61 = 3.0835 
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Figure 4: The storage modulus E' MPa versus the frequency of oscillations u Hz for 
carbon black filled rubber in tensile dynamic tests on preloaded specimens (e = —0.1) 
with the amplitude of oscillations Ae = 0.006 at temperature T K. Circles: experimental 
data p3| . Solids lines: results of numerical simulation. Curve 1: T — 253; curve 2: 
T = 296; curve 3: T = 333; curve 4: T = 373 
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Figure 5: The rate of breakage T s 1 versus the degree of undercooling AT K for a carbon 
black filled rubber. Circles: treatment of observations 
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Solid lines: approximation of 
the experimental data by Eq. (64) with do = —2.5489 and d\ = 0.0272 
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